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ABSTRACT 


BIGRAP  is  a  package  of  interactive  graphics  programs, 
written  for  use  on  the  TEKTRONIX  4014  Graphics  Terminal 
connected  to  the  ARRADCOM  CDC  6500/6600  computer  configura¬ 
tion.  This  package  consists  of  a  set  of  intricate  programs 
that  allows  a  user  to  input  component  success/failure  data 
and  a  Boolean  expression  depicting  system  reliability  logic  for 
the  purpose  of  assessing  system  reliability.  The  computer 
converts  the  logic  expression  to  an  algebraic  expression  for 
the  system  reliability  as  a  function  of  the  individual  com¬ 
ponent  reliabilities.  A  Bayesian  statistical  algorithm  is 
then  employed  to  provide  the  user  with  point  and  confidence 
interval  estimates  of  system  reliability.  In  addition  the 
graphics  feature  of  the  package  displays  histograms  and 
corresponding  Beta  distributions  involved  in  the  analysis. 

The  programs  are  run  sequentially  as  a  cataloged  procedure 
using  the  BEGIN/REVERT  utility.  This  report  is  intended  to 
serve  as  a  guide  to  potential  users  of  BIGRAP  for  operation 
and  interpretation  of  results. 
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2.  INTRODUCTION  AND  METHODOLOGY 


The  Bayesian  Interactive  Graphics  Reliability  Assessment 
Procedure  (BIGRAP)  is  a  package  of  interactive  graphics  programs 
written  for  use  on  the  Tektronix  4014  Graphics  Terminal  tied  into 
the  ARRADCOM  CDC  6500/6600  computer  system.  The  package  consists 
of  a  set  of  Fortran  programs  which  are  executed  sequentially  as 
a  catalogued  procedure  using  the  BEGIN/REVERT  utility.  Each  of 
the  five  main  programs  in  BIGRAP,  along  with  its  associated 
subroutines,  performs  a  major  computational  step  required  as 
part  of  a  pseudo-Bayesian  approach  to  reliability  assessment  of 
complex  systems  based  on  component  success/failure  data.  This 
report  is  intended  to  serve  as  a  BIGRAP  user  manual,  as  well  as 
an  introductory  methodological  summary  providing  the  reader  with 
the  essentials  necessary  to  use  BIGRAP  and  understand  the  output 
it  generates. 

Reliability  assessment  approaches  for  complex  systems 
have  historically  had  two  major  weaknesses.  Firstly,  they 
provided  only  a  point  estimate  of  reliability  and  failed  to 
provide  a  measure  of  the  total  uncertainty  reflected  by  the 
limited  component  data.  Secondly,  the  computational  require¬ 
ments  were  unwieldy  and  oriented  to  individual  systems 
requiring  a  large  amount  of  time  to  accomplish  with  minimal 
computer  aid.  BIGRAP  has  neither  of  these  weaknesses.  A 
Bayesian  reliability  assessment  approach  is  employed  for  com¬ 
puting  point  estimates,  as  well  as  confidence  limits  providing 
the  required  measures  of  uncertainty.  The  interactive 
graphics  application  and  BEGIN/REVERT  utility  make  the 
procedure  easy  to  use,  even  by  users  unfamiliar  with  computer 
programming. 

In  this  section,  the  approach  underlying  BIGRAP  will  be 
outlined  so  that  the  potential  user  will  have  some  idea  of 
(1)  the  preliminary  work  he  will  have  to  do  before  running 
BIGRAP  to  set  up  his  problem;  (2)  *  the  rationale  behind  the 
computations  performed  by  BIGRAP;  and  (3)  the  meaning  of 
the  output  generated  by  BIGRAP. 

The  first  step  in  the  system  reliability  assessment  approach 
underlying  BIGRAP  is  the  determination  of  the  total  number  of 
components  and  component  types  which  make  up  the  system. 

Components  are  considered  to  be  of  the  same  type  in  this  context 
if  they  are  identical  in  design  and  manufacture,  i.e.  they  can 
be  considered,  in  the  statistical  sense,  to  be  drawn  from  the 
same  population.  The  meaning  and  importance  of  the  distinction 
between  the  total  number  of  components  and  the  number  of  component 
types  will  become  clear  to  the  reader  as  he  proceeds. 
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The  next  step  in  the  assessment  procedure  is  the  construction 
of  a  reliability  block  diagram  which  depicts  the  logic  by  which 
system  success  (or  failure)  occurs  as  a  result  of  success 
(or  failure)  combinations  of  all  the  components  which  comprise 
the  system.  While  it  is  likely  that  most  potential  users  of 
BIGRAP  are  familiar  with  the  concept  of  a  reliability  block 
diagram,  a  simple  example  will  be  provided  here  so  that  succeeding 
steps  of  the  assessment  approach  may  be  more  easily  illustrated. 
Consider  the  simple  system  functional  schematic  provided  in 
Figure  1. 


3  MOVEMENT 
MECHANICAL  TIMER 


Figure  1 
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RELIABILITY  BLOCK  DIAGRAM 


Component  Name 


Component  Component 
Type 

Symbol_ Symbol 


Component 

Type 

Number 


Acceleration  Switch 

A 

A 

1 

Battery  1 

Bl 

B 

2 

Battery  2 

B2 

B 

2 

Timer  Movement  1 

T1 

T 

3 

Timer  Movement  2 

T2 

T 

3 

Timer  Movement  3 

T3 

T 

3 

Detonator 

D 

D 

4 

Explosive  Charge 

E 

E 

5 

Figure  2 
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This  simple  system  represents  the  author's  naive  conception 
of  a  possible  (but  improbable)  missile  high  explosive  warhead 
section.  The  system  is  designed  so  that  upon  sensing  sufficient 
acceleration  after  launch  the  acceleration  switch  closes  and  the 
three  identical  movements  of  the  timer  are  started.  When  a 
preset  time  is  measured  by  any  one  of  the  three  timers  the 
circuit  is  closed  between  the  detonator  and  the  two  identical 
batteries,  each  having  sufficient  voltage  to  fire  the  detonator. 
The  firing  of  the  detonator  ignites  the  explosive  charge 
completing  the  system  mission.  The  system  has  a  total  of  eight 
components,  but  only  five  component  types.  Figure  2  lists  the 
different  components  along  with  their  assigned  component  symbols 
and  type  symbols  below  the  reliability  logic  block  diagram  for 
the  system.  The  block  diagram  reflects  that  components  A,  D,  E, 
one  or  both  of  components  Bl  and  B2 ,  and  at  least  one  of 
components  Tl,  T2  and  T3  must  function  successfully  to  achieve 
system  success. 

BIGRAP  was  developed  with  the  intention  of  allowing  the 
user  to  provide,  as  part  of  the  input,  the  system  logic  as 
depicted  by  the  reliability  block  diagram.  Boolean  algebra 
provides  the  tools  with  which  this  can  be  accomplished.  Once 
the  reliability  block  diagram  has  been  constructed  it  is  easy 
to  formulate  the  equivalent  Boolean  algebraic  expression.  For 
the  sample  system  of  Figure  1  the  Boolean  expression  for  system 
success,  S,  is 


S=A* (Bl  +  B2)*((T1  +  T2 )  +  T3)*D*E 


where  the  symbols  "*"  denotes  "and"  logic  and  "+"  denotes  "or" 
logic.  If,  in  this  expression,  one  assigns  to  the  component 
symbols  1  for  success  or  0  for  failure  and  then  evaluates  the 
expression  according  to  the  following  simple  rules: 


0*0  =. 0 
1*1  =  1 


0*1  =  1*0  =  0 
1+1  =  1 
0+1  =1+0=1 
0+0  =  0 


(Note:  The  usual  associative, 

distributive  and  commutative 
axioms  hold  as  well) 


a  value  of  S=1  or  S=0  will  result  in  implying  that  system  success 
or  failure,  respectively,  is  achieved  by  the  assigned  component 
successes  and/or  failures.  For  example,  if  one  assigns  A=l, 

Bl=0,  B2=l,  Tl=l,  T2=l ,  T3=0,  D=1 ,  and  E=1  the  Boolean  expression 
with  these  values  substituted  becomes 


S  =  1*(0+1)*( (1+1)+0)*1*1 
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which  can  be  reduced  via  the  above  rules  to 


S  =  1*1* (1+0) *1*1 


and  further  to 


S  =  1*1*1*1*1  =  1 

representing  system  success.  As  another  example  consider  the 
assignment  A=l,  Bl=0 ,  B2=0,  Tl=l,  T2=0,  T3=0,  D=l,  E=1  resulting 
in 


S  =  l*(0+0)*((l+0)+0)^l*l 
=  1*0*(1+0)*1*1 
=  1*0*1*1*1 
=  0 

showing  that  this  latter  assignment  results  in  system  failure. 

For  readers  for  which  the  above  explanation  does  not  sufficiently 
clarify  the  process  of  translating  the  block  diagram  to  a  Boolean 
expression,  many  papers  and  texts  in  the  extant  literature  provide 
more  detailed  explanations  (see  references  11  and  12). 

Boolean  algebra  provides  the  means  by  which  the  system  event 
can  be  described  in  terms  of  component  events.  The  subject  of 
reliability,  however,  deals  with  probabilities,  and  in  particular 
the  probability  of  the  system  success  event.  In  fact  the  system 
reliability  is  defined  here  as  the  probability  of  system  success. 
System  reliability  assessment  based  on  component  data  involves 
the  computation  of  point  and  confidence  interval  estimates  of 
system  reliability  using  success/failure  data  which  directly 
estimates  the  individual  component  reliabilities.  The  component 
reliabilities  are  defined  analogously  as  the  probabilities  of 
component  success.  Components  of  the  same  type  have,  by 
definition,  the  same  reliability  value.  To  proceed  with  the 
reliability  assessment  approach,  therefore,  a  procedure  is 
required  by  which  the  Boolean  logic  expression  can  be  converted 
to  an  algebraic  expression  for  system  reliability  in  terms  of 
the  component  type  reliabilities.  The  first  three  programs 
executed  as  part  of  BIGRAP  perform  such  a  conversion  for  a 
Boolean  expression  of  a  system  provided  as  input.  The  method 
employed  therein  to  accomplish  this  is  too  lengthy  to  describe 
here  but  is  adequately  explained  in  References  n  and  12. 

The  user  need  only  be  acquainted  with  the  form  of  the  results 
which  will  be  illustrated  herein. 
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The  conversion  is  accomplished  in  two  stages.  The  first 
stage  involves  determination  of  the  total  number  of  success 
paths  (or  tie  sets).  For  our  example  there  are  six  such 
success  paths.  They  are: 

A*B1*T1*D*E 

A*B1*T2*D*E 

A*B1*T3*D*E 

A*B2*T1*D*E 

A*B2*T2*D*E 

A*B2*T3*D*E 

The  second  stage  uses  this  essentially  expanded  form  of  the 
Boolean  expression  and  converts  it  to  an  ordinary  algebraic 
expression  for  the  system  reliability,  Rg;  in  terms  of  the 

component  type  reliabilities,  resulting  in 

RS  -  6rarbrtrdre  -  6rarb4rdre 
+  2rarbrtrdre  -  3rarbWe 

+  3rarXrdre  -  rar^rdre 

where  an  on  the  right  side  of  the  equal  sign  denotes  the 

reliability  of  component  type  X.  This  expression  provides  the 
required  relationship  between  the  system  reliability  and  the 
component  type  reliabilities. 

The  Bayesian  aspects  of  the  assessment  procedure  are 
introduced  at  this  point.  The  basic  underlying  assumption  is 
that  the  unknown  component  type  reliabilities  (measured  by  the 
success/failure  data)  can  be  considered  as  independent  random 
variables  in  the  Bayesian  sense,  each  having  a  prior  probability 
distribution  belonging  to  the  family  of  beta  distributions. 

The  general  form  of  the  beta  probability  density  function  (pdf) 
for  a  random  variable,  R,  is: 

f  g  (R  | A, B)  =  r  (AM-CB)  RA"1(l-R)B"1for  0  <  R  <  1 


=  0  otherwise. 

where:  (i)A  and  B  are  the  parameters  of  the  beta  distribution 

(with  the  restriction  that  A  >  0  and  B  >0. 
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(ii)T(.)  is  the  complete  gamma  function  defined  by 


r  (x)  =i[yX'1e"ydy 


having  the  property  that  r (x+1)  =  xr(x)  so  that  for  a  positive 
integer,  n, 

r(n+l)  =  n! 

If  the  component  types  are  indexed  1  thru  m,  where  m  is  the 
number  of  types  (m  =  5  for  the  example)  the  prior  beta  pdf  of 
each  component  type  reliability  will  be  represented  by 
fg(Ril  ,  i=l,...,m.  BIGRAP  allows  the  user  to  select  from 

a  variety  of  options  for  defining  the  A^  and  B^  values.  These 

options  will  be  described  in  Section  5b  of  this  report. 

The  data  is  introduced  into  the  assessment  procedure  by 
similarly  indexing  the  component  type  number  of  successes 
and  number  of  failures  as  s^  and  f ^  respectively.  For  a  given 

component  type  reliabilities,  Ri ,  i=l,...,  m,  the  likelihood  of 

observing  the  given  data  s^,  f^  is  then  provided  by  the  binomial 

distribution  with  probability  mass  function  given  by 


The  likelihood  expressed  as  a  function  of  R^  is  then 


Application  of  Bayes  Theorem  (Principle  of  Inverse  Probability) , 
which  says  that  the  posterior  (updated  with  data)  probability  is 
proportional  to  the  product  of  the  prior  probability  and  the 
likelihood,  results  in  the  following: 


Posterior  pdf  of  R^ 


0 


The  expression  on  the  riaht  can,  upon  substitution  of  the  com¬ 
ponent  type  prior  beta  pdf  and  binomial  likelihood,  be  reduced  to 
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so  that  the  posterior  distribution  of  Ri  is  also  a  member  of  the 
beta  family  with  posterior  parameters  given  by 


AX.  =  A.  +  s. 
ill 


BX.  =  B.  +  f. 
111 


It  is,  therefore,  a  simple  matter  to  add  the  component  success/ 
failure  data  into  the  assessment  procedure  for  a  given  set  of 
component  type  prior  distributions. 

Since  each  of  the  component  type  reliabilities  is  considered 
a  random  variable,  the  system  reliability,  Rs,  which  is  a  function 

of  the  component  type  reliabilities,  is  also  a  random  variable. 

The  exact  form  of  the  probability  distribution  of  R  clearly 

depends  on  the  distribution  of  the  component  type  reliabilities 
and  the  relationship  of  the  component  type  reliabilities  to 
system  reliability  as  reflected  in  the  algebraic  reliability 
expression  for  Rg.  In-depth  investigations  have  shown  that 

this  exact  distribution  can  be  reasonably  approximated  by  a 
member  of  the  beta  family. 

BIGRAP  uses  two  methods  for  determining  the  appropriate 
beta  fit  to  the  distribution  of  Rs.  Monte  Carlo  sampling  and 

the  method  of  moments.  The  Monte  Carlo  approach  involves 
generation  of  a  large  number  of  "trials",  where  on  each  trial 
each  of  the  component  type  reliability  posterior  distributions 
f  g(Ril  AXjyBX^),  i  =  l,...,m,  is  randomly  sampled  yielding  a 

set  of  type  reliability  values  R.,R,,...,R  .  These  values  are 

la  m 

substituted  into  the  algebraic  reliability  expression  and  a 
value  of  Rs,  the  system  reliability,  is  obtained.  This  process 

is  repeated  for  the  total  number  of  trials  specified  by  the  user. 

A  histogram  depicting  the  distribution  of  the  resulting  system 
reliability  values  is  constructed.  In  addition  the  arithmetic 
—  f  o 

mean,  Rg,  and  sample  variance,  s^,  of  the  Rg  values  is  computed. 

A  direct  relationship  between  the  parameters  of  a  beta  distribu¬ 
tion  and  its. mean  and  variance  is  then  used,  with  R  substituted 

s 

2 

for  the  mean  and  s  substituted  for  the  variance,  to  estimate 

K 

the  parameters  of  a  system  beta  distribution  serving  as  a  beta 
fit  to  the  histogram. 

(  In  discussing  the  method  of  moments  a  few  basics  should 

be  recognized  concerning  the  moments  of  continuous  random 
variables  in  general  and  beta  distributed  random  variables  in 
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particular.  It  is  assumed  that  the  reader  is  familiar  with 
the  concept  of  expectation.  The  rth  moment  of  a  random 
variable,  X,  is  defined  to  be  the  expectation  of  the  rth 
power  of  X,  i.e. 

th  t* 

r  moment  of  X  =  E(X  ) 

The  first  moment  is  referred  to  as  the  mean,  y,  of  the 
•  .  •  0 
distribution  of  the  random  variable.  Its  variance,  a  ,  is 

defined  as  the  second  moment  about  the  mean.  Therefore, 

H  =  E  (X) 
a2  =  E  ( X-  v  )  2 
or  equivalently 

ct2  =  E(X2)-E(X)2 

This  latter  form  of  the  expression  for  o2  follows  from  the 
following.  If  X  and  Y  are  random  variables  (or  functions  of 
random  variables)  and  and  C constants,  then 

e(c1x+c2y)  =  c1e(x)+c2e(y) 


If  X  and  Y  are  independent,  then 
E  (XY)  =  E  (X)  E  (Y) 


For  a  random  variable  R  having  a  beta  distribution  with  pdf 
given  by  F  ^(R|  A,B)  the  r*"*1  moment  is 


E(Rr) 


r  A+i-1 


A+B+i-1 

i=l 


This  result  then  yields  the  following  expression  for  m  and  a2 
of  the  beta  distribution 


M  =  A/ (A+B) 

2  _  A  (•  A+l  A  ^ 

°  A+B  *>  A+B+l  A+B  > 


12 


By  solving  these  two  expressions  for  A  and  B  in  terms  of  v 
2 

and  cr  the  following  useful  formulae  are  obtained: 


A  = 


B  = 


2(1-  U  )  _ 

„  2 


p 


1 


2  . 

By  substituting  estimates  of  p  and  a  into  these  last  two 
expressions,  estimates  of  A  and  B  can  be  obtained.  This  is 
precisely  what  was  done  in  determining  the  beta  fit  to  the 
.  .  —  2 

Monte  Carlo  histogram  wherein  R  and  s_  were  used  as  estimates 
2 

of  n  and  a,  respectively. 

The  method  of  moments  involves  the  computation  of  the 
first  and  second  moments  of  the  posterior  distribution  of 
system  reliability  in  terms  of  the  moments  of  the  component 
type  reliability  posterior  distributions.  The  algebraic 
reliability  expression  for  R  and  its  algebraic  square  for 
2  s 

Rg  are  used  for  this  purpose.  The  example  system,  with  the 
type  numbers  listed  in  Figure  2,  can  be  used  to  illustrate 
this  procedure.  The  first  moment  of  R  is  first  determined 
as :  s 


E  ( R  )  =  E  ( 6R.  R„R-R .  Rc  —  6R.R0R-  R.R,, 
s  12345  12345 


+2R1R2R3j5R4R5  -  etc.) 


Using  what  has  been  stated  about  expectations,  this  can  then 
be  expressed  as 

E  (Rg)  =  6E(R1)E(R2)E(R3)E(R4)E(R5) 

-  6E(R1)E(R2)E(R32)E(R4  )E(R5) 

+  2e(r1)e(r2)e(r33)e(r4)e(r5) 


-  etc. 

i 

By  substituting  the  appropriate  expression  for  each  beta 
moment  in  terms  of  the  component  type  posterior  beta 
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The  appropriate 


parameters,  E(R  )  can  be  evaluated. 

o 

substitutions  follow  from; 


r  r  AX.+j-l 

E (Ri  *  “  j=l  AX.+BX.+j-l 


l  l 
AX, 


e.g.  E (R^)  AX,+BX, 


,  e(r3  ) 


3 

E(R3  >  =  AXT+fc: 


ax3+i 


AX3+BX3+1 


AX3  ax3+1 

AX3+BX3  *  AX3+BX3+1 

AX3+2 

AX3+BX3+2 


.  .  2 

Similarly  E(Rg  )  is  evaluated  using  the  expectations  of  the 

algebraic  expression  for  the  square  of  R  ,  i.e., 

s 


E(R  2)  =  E(36R  2R  2R  2R  2R  2  -  72R  2R  2R  3R  2R  2 
s  12345  12345 


+  24R2R  2R  4R  2R  2  -  36R.2R  3R  2R  2R  2 
12345  12345 

+  etc.) 


By  setting 


U  =  E  ( R  ) 
s  s 


and 


=  E(R  2)  -  E  ( R  )  2 

o  S 


in  the  expressions  for  the  beta  A  and  B  parameters,  the  method 
of  moments  posterior  beta  distribution  for  system  reliability, 
Rsf  is  fully  determined  by  the  resultant  A  and  B  values;  i.e. 


A 


s 


B  = 
s 


+  ys“1 


f 


The  result  obtained  by  either  the  Monte  Carlo  method  or 
the  method  of  moments  is  a  beta  distribution  which  approximates 
very  closely  the  true  posterior  distribution  of  the  system 

reliability,  Rg.  In  each  case  the  resultant  beta  distribution 
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is  completely  defined  by  its  two  parameters,  A  and  B,  which 

be  denoted  by  A^  and  Bg  since  they  are  the  parameters  of 

the  posterior  distribution  of  R  .  The  user,  however,  is 

s  ’ 

generally  performing  the  assessment  procedure  for  the  purpose 
of  obtaining  point  or  confidence  interval  estimates  of  the 
system  reliability.  BIGRAP  extracts  these  estimates  from  the 
information  contained  in  the  system  posterior  beta  distribution. 
In  all  cases,  the  mean  and  mode  (if  it  exists  and  is  between  0 
and  1.0)  of  the  posterior  distribution  of  R  are  provided.  These 
values  are  computed  as  follows:  s 

Mean  of  R  =  A  /(A  +B  ) 
s  s  s  s 

If  A_  >  1  and  B  >  1,  then 
s  s 

Mode  of  Rs  =  ( Ag-1) / (Ag+Bg-2 ) 

The  recommended  strategy  for  point  estimation  is  to  take  as  the 

point  estimate  of  R  tne  mode  when  A  and  B  are  both  greater 

_ s _  s  s 

than  (not  equal  to)  1.0  and  the  mean,  otherwise" 


For  the  purpose  of  confidence  interval  estimation,  BIGRAP 
presently  provides  the  lower  50,  90,  and  95  per  cent  Bayesian 
confidence  limits.  The  lower  100  (l-y)%  confidence  limit  for 
Y=  .5,  .1,  or  .05  is  defined  as  the  value  of  R  satisfying  the 
following: 


where  the  A  and  B  that  are  used  result  from  either  the  Monte 

b  S 

Carlo  method  or  the  method  of  moments.  Precise  numerical 
integration  routines  included  in  BIGRAP  set  up  and  interpolate 
a  table  of  y  vs  R  to  determine  the  appropriate  value  of  R. 

For  a  particular  system  configuration  to  which  BIGRAP  is 
being  applied  for  the  first  time  the  user  may  wish  to  employ 
both  the  Monte  Carlo  approach  and  the  method  of  moments.  By 
comparing  the  point  and  confidence  limit  estimates  resulting 
from  the  both  methods  the  user  can  satisfy  himself  that  the 
assumptions  involved  in  the  method  of  moments  approach  (use  of 
only  first  two  moments  and  a  beta  fit  to  system  reliability 
distribution)  do  not  introduce  significant  errors.  In  subsequent 
applications  to  the  same  system  configuration  the  user  need  only 
apply  the  method  of  moments  should  the  comparison  be  favorable. 
If,  on  the  other  hand,  there  is  significant  disagreement  between 
these  results,  or  if  the  beta  distribution  in  the  Monte  Carlo 
approach  does  not  appear  to  fit  the  histogram,  the  user  may 
then  choose  to  continue  with  the  Monte  Carlo  approach.  In 
the  latter  case  the  user  should  make  use  of  the  histogram 
itself,  rather  then  its  beta  fit,  to  extract  empirical  point  and 
interval  estimates. 
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The  succeeding  sections  of  this  report  present  a  detailed 
description  of  the  operation  and  structure  of  the  BIGRAP 
software.  Instructions  for  job  control  and  input  are  also 
provided.  Several  examples  of  the  use  of  BIGRAP  are  presented. 

In  addition  the  complex  manual  computations,  which  are  the 
alternative  to  employing  BIGRAP,  for  a  simple  system  configuration 
are  exhibited  to  emphasize  the  value  of  the  BIGRAP  package. 

The  reader  interested  in  more  details  regarding  the 
theoretical  basis  for  BIGRAP  should  look  to  the  cited  references 
and  in  particular  Reference  14. 

As  a  final  introductory  comment,  the  authors  wish  to  point 
out  that  a  large  amount  of  work  was  performed  in  comparing  the 
results  from  BIGRAP  to  other  approaches  to  reliability  assessment 
both  classical  and  Bayesian.  For  the  restricted  cases  (in  the 
sense  of  system  configuration  complexity)  to  which  these  other 
approaches  could  be  applied,  BIGRAP  performed  at  least  as 
well  and,  in  most  cases,  better.  So  whether  the  reader  considers 
himself  (or  herself)  a  Bayesian,  Classicist,  or  "middle-of-the- 
roader"  the  assessment  algorithm  of  BIGRAP  should  prove 
satisfactory. 


16 


3.  DESCRIPTION  OF  THE  MAIN  PROGRAMS 


The  Bayesian  Interactive  Graphics  Reliability  Assessment 
Procedure  is  made  up  of  five  main  computer  programs  (B0LN1 , 

B0LN2 ,  B0LN3 ,  MOMENT,  MONTI)  which  will  be  described  here: 

Program  B0LN1  -  Uses  as  input,  the  Boolean  expression 
of  a  reliability  block  diagram  in  free  form.  The  program 
transforms  the  simple  expression  to  integer  form  by  comparing 
each  element  in  the  expression  to  the  list  of  special  characters 
(plus  sign,  asterisk,  and  right  and  left  parentheses)  and 
components  and  places  the  number  of  the  list  that  equals  the 
element  of  the  expression  in  the  same  position  as  the  element. 

The  program  then  goes  through  the  expression  and  stores  the  number 
of  sections  and  the  locations  of  the  plus  signs  and  left 
and  right  parentheses  for  each  set  of  parentheses.  The  results 
of  the  above  operations  are  written  on  a  disk  for  use  in  program 
B0LN2 . 

Program  B0LN2  -  Reads  the  disk  output  from  program  B0LN1 
as  input.  The  program  expands  the  expression  in  the  conventional 
manner  to  ascertain  all  the  success  paths.  As  each  path  is 
derived,  it  is  compared  to  all  previously  generated  paths.  If 
the  newly  derived  path  contains  or  is  equal  to  any  previously 
generated  path,  the  new  one  is  eliminated  as  a  redundancy.  If 
the  newly  derived  path  is  contained  in  any  other  previously 
generated  path,  the  new  one  replaces  the  older  one,  which  is 
then  eliminated  as  a  redundancy.  If  neither  of  the  above  is  found 
to  be  true,  the  newly  ascertained  path  is  added  to  the  list  of 
non-redundant  paths.  The  results  are  then  written  (on  both  disc 
and  printout) . 

Program  B0LN3  -  Converts  Boolean  logic  expression  to  an 
expanded  algebraic  reliability  expression.  Constructs  a  Fortran 
subroutine,  on  TAPE  12,  for  calculation  of  system  reliability 
using  the  algebraic  expresison,  passes  this  subroutine  to  program 
MONTI,  passes  the  algebraic  reliability  expression  to  output  and 
program  MOMENT  on  TAPE  7 . 

Program  MOMENT  -  Allows  the  user  to  select  one  of  four 
options.  Option  1  uses  the  algebraic  reliability  expression 
from  program  B0LN3.  It  then  determines  and  stores  the  square 
of  the  algebraic  expression  and  formulates  the  system  distri¬ 
bution  1st  and  2nd  moments  as  functions  of  the  moments  of  the 
component  distributions.  These  results  are  then  used  to  determine 
the  equal  component  prior  beta  parameters  A'  and  B' ,  which  yield 
a  uniform  prior  distribution  for  the  system  reliability.  It  also 
computes  the  posterior  beta  parameters  for  each  component  and 
substitutes  these  results  into  the  expressions  for  the  system 
moments.  The  method  of  moments  estimates  of  system  posterior 
beta  parameters  are  then  computed. 
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In  Option  2  and  Option  4  the  user  provides  the  prior  beta 
parameters  (either  equal  or  unequal)  for  each  component.  There 
is,  of  course,  no  need  for  a  solution  to  the  moment  equations 
prior  to  performing  the  method  of  moments  estimation  procedure. 

In  Option  3  the  moment  expressions  are  set  equal  to  user 
specified  1st  and  2nd  prior  moments.  These  equations  are  then 
solved  and  the  same  procedures  are  followed  as  in  Option  1. 

Program  MONTI  -  Uses  the  component  posterior  beta  parameters 
for  each  component  as  well  as  the  method  of  moments  estimates 
of  the  system  beta  parameters  as  input.  The  component  parameters 
are  used  to  define  beta  distributions  for  each  component.  Monte 
Carlo  simulation  is  then  performed  by  sampling  in  turn  each  of 
the  component  distributions  for  a  component  reliability  value; 
substituting  these  values  into  the  system  algebraic  reliability 
expression  (generated  in  B0LN3) ,  generating  a  system  reliability 
value  on  each  Monte  Carlo  trial.  After  a  sufficiently  large 
number  of  trials,  a  histogram  of  system  reliability  values  is 
obtained,  in  addition  to  the  mean  and  variance  of  the  generated 
system  reliability  values.  It  then  computes  lower  confidence 
limits  for  system  reliability  based  on  the  method  of  moments  and 
Monte  Carlo  results.  Finally  the  histogram  is  displayed  on  the 
Tektronix  screen  with  a  beta  curve  fit  from  the  mean  and  variance 
superimposed.  A  useful  feature  of  program  MONTI  is  the  capability 
of  generating  additional  histograms  for  the  same  set  of  inputs 
varying  either  the  number  of  histogram  cells  (increments)  or  the 
number  of  Monte  Carlo  trials  (samples) .  The  random  number 
generator  subroutines,  RDMIN  and  RDMOUT,  are  also  employed  to 
avoid  the  use  of  the  same  or  overlapping  set  of  random  numbers 
on  subsequent  runs  with  the  same  set  of  inputs.  These  features 
are  illustrated  in  section  5  of  this  report. 

The  present  form  of  Programs  B0LN1 ,  B0LN2,  and  B0LN3  were 
based  on  programs  originally  written  by  Howard  Corneilson.  Each 
of  his  original  separate  programs  were  run  as  batch  jobs.  The 
modified  programs  are  used  interactively  with  the  graphics  terminal 
and  are  run  sequentially  using  BEGIN/REVERT. 

Program  MOMENT  was  based  on  a  batch  program  written  by 
Bruce  D.  Barnett.  The  BIGRAP  form  of  the  program  derives  its 
input  directly  from  program  B0LN3  without  the  need  for  the  user 
to  provide  the  algebraic  expression. 

Program  MONTI  was  based  on  two  programs  written  by  Irving 
E.  Teviovitz  for  the  274/1700/6500  large  scale  refresh  graphics 
facility.  One  program  illustrated  a  user  specified  beta  distri¬ 
bution  probability  density  function  and  the  second  program  il¬ 
lustrated  a  beta  random  distributed  histogram.  The  two  programs 
were  combined  and  converted  for  use  on  the  Tektronix  terminal.  The 
JBETA  subroutine  in  MONTI  for  numerical  integration  of  the  beta 
distributions  was  developed  by  Joseph  E.  Bay  through  modification 
of  subroutines  used  in  program  RELY  (Neff  1969) . 
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The  BEGIN/REVERT  utility,  which  automatically  executes 
all  the  job  control  statements  that  are  needed  to  run  these 
programs,  is  a  set  of  programs  developed  by  the  University 
of  Washington. 
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4.  TERMINAL  INITIALIZATION  AND  LOGIN 


A  typical  Tektronix  4014  Graphics  Terminal  set-up  is 
pictured  in  Fiaure  3.  A  complete  description  of  the  terminal 
and  its  operation  can  be  obtained  in  Reference  3. 

The  purpose  of  this  section  is  to  accruaint  the  potential 
user  with  some  basic  information  required  to  operate  the  4014 
terminal  in  connection  with  use  of  BIGRAP. 

Figure  4  orovides  a  description  of  the  keyboard  con¬ 
figuration  of  the  4014  terminal. 
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Typical  Tektronix  4014 

Graphics  Terminal 
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This  is  the  Keyboard  of  the 
4015  Computer  Graphics  Terminal 


The  following  procedures  are  required  for  Terminal 
initialization. 

1.  Place  the  LOCAL/LINE  switch  to  LOCAL  (the  terminal  will 

be  isolated  from  the  computer  and  causes  keyboard  inputs  to 
be  executed  by  the  terminal) . 

2.  Turn  the  Terminal  ON  (the  power  switch  is  on  the  front 
lower  right  of  the  console  (one  foot  from  floor) .  The 
green  power  light  on  the  left  of  the  keyboard  will  light 
and  the  display  screen  will  become  bright  after  a  few 
seconds . 

NOTE:  At  turn  on,  the  terminal  is  in  Alpha  Mode  with  the 

cursor  at  HOME  position  with  normal  (large  character 
size  and  is  in  synchronous  communications  mode. 

3.  Turn  on  hard  copier  (if  one  is  available).  Erase  one 
display  by  pushing  the  RESET-PAGE  key  (at  least  one  erase 
after  power-up  is  required  to  start  the  interface) . 

4.  Set  the  following  keyboard  switches: 

CODE  EXPANDER  TO  OFF 
CLEAR  WRITE  TO  OFF 

5.  Set  the  interface  control  switches 

Rotary  baud  switch  to  300  for  300  baud  rate  (rear  of 
terminal)  ASCII/BCD  switch  to  ASCII  (rear  of  terminal) 

6.  SHIFT,  RESET-PAGE  keys  are  depressed  on  the  keyboard 
simultaneously.  This  will  initialize  the  Terminal,  by 
placing  initial  conditions  on  the  Terminal  circuits. 

7.  Press  CTRL,  SHIFT  and  P  keys  -  (for  local  echo). 

8.  Press  the  following  three  keys  sequentially: 

RETURN  ESC  :  (For  small  characters) 

9.  Put  LOCAL/LINE  switch  to  LINE 

10.  Remove  grey  phone  from  the  modern  cradle,  switch  phone  to 
TALK.  Phone  select  switch  under  table  to  be  in  clockwise 
position  marked  "OPT’’. 

11.  Dial  up  (to  connect  to  CDC  6500/6600) 

12.  When  the  computer  gives  an  audible  tone,  the  button 
marked  DATA  is  pushed  and  the  phone  is  hung  up. 
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13.  INTERCOM  sign-on  (LOGIN)  procedures  will  then  be  re¬ 
quested  on  the  terminal  screen  and  be  responded  to  by 
the  user  as  illustrated  in  figure  5. 
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5.  ILLUSTRATIVE  EXAMPLE 


INTERCOM  commands  to  initiate  the  B10RAP  procedure  also  il¬ 
lustrated  in  the  same  fiqure. 

To  show  what  happens  next,  the  example  from  the  Introduction 
(Section  2)  will  be  used. 

The  screen  is  cleared  and  information  shown  on  page  27  is 
requested.  The  user's  inputs  are  underlined. 

The  name,  type,  number  of  successes  and  failures  of  each 
component  is  entered  one  at  a  time,  only  when  prompted  and 
each  followed  by  a  return. 

Note  the  warninq  for  entering  the  Boolean  expression. 

At  the  end  of  the  page  a  Copy  may  be  made  and/or  the  "ESC" 
key  followed  by  the  "Return"  key  must  be  depressed  (this 
sequence  must  be  followed  at  the  end  of  every  page) . 

Systems  messages  are  printed  out  on  page  30. 

The  output  from  BOLN  2  and  BOLN  3  are  batched  to  terminal 
specified  above. 

This  output  consists  of  system  success  paths  shown  on  page 
28  and  algebraic  reliability  expression,  pace  29. 

Various  options  are  requested  on  page  31. 

A  "Copv , Esc , Return- (C ,E , R) "  sequence  is  required  here.  Various 
results  are  shown  on  page  32. 

A  "C,E,R"  sequence  is  used  at  end  of  "THE  VAR  OF  RS=  ". 

On  page  33  RANDOM  START  is  any  ten  diait  inteaer  number.  The 
NUMBER  of  Increments  and  Samples  are  floating  point  numbers. 

A  "C,E,R,"  sequence  is  used  at  end  of  "RANDOM  END  IS  ". 

The  histogram  and  distribution  curve  is  shown  on  page  34. 

A  "C,E,R"  is  required  at  the  end  of  cur  e.  Next  follow 
directions  on  oaae  34. 

Pages  3  5  to  38  show  how  the  histogram  changes  as  the  number 
of  increments  change  (line  2  of  page  35  and  37.) 

The  display  on  page  35  is  shown  when  the  cursor  is  placed  on 
line  3  and  the  number  "3"  is  pressed. 

Panes  39  to  50  show  how  Options  2,3  and  4  are  used. 
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CONTROL  DATA  INTERCOM  4.5 
DATE  06/08/78 
TIME  09.56.41. 
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Figure  5 

LOGIN  Procedures  (User  Inputs  are  Underlined) 
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place  cursor  on  this  line  and  press  proper  KEV  for  desired  action 

ICEV  ACTION 

1  JO»  TERN I NAT I ON 

2  NEXT  CASE-NEU  A  AND  B 

3  NEXT  CASE-OLD  A  AND  I 
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PLACE  CURSOR  ON  THIS  LINE  AND  PRESS  PROPER  KEY  FOR  DESIRED  ACTION 
KEY  ACTION 

1  JOB  TERMINATION 

2  NEXT  CASE-NEU  A  AND  B 

3  NEXT  CASE-OLD  A  AND  B 
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PLACE  CURSOR  ON  THIS  LINE  AND  PRESS  PROPER  KEY  FOR  DESIRED  ACTION 
KEY  ACTION 

1  JOB  TERMINATION 

2  NEXT  CASE-NEU  A  AND  B 

3  NEXT  CASE-OLD  A  AND  B 


38 


4-> 

C 

Q) 


0 

to 

i — 

z 

K 

Ui 

l  , 

e 

*■1 

>T. 

o 

0 

V 

to 

a, 

eq; 

C 

tOUJUi 

rH 

0:1-1- 

UitOUl 

1 —  >E 

c 

UitO<E 

0 

e  a: 

*  i 

4-> 

<ro<r 

a, 

CLZCL 

o 

<r  o 

hoq: 

'h 

j 

UIO 

0 

a:to»-« 

Cl4 

o  on 

►— « Q  Q_ 

o:z 

0 

cl  <r  —> 

CD 

<r 

tn 

CE—JI-Z) 

D 

o<r  too 

*-hZ)CZUJ 

aOMZ 

J3 

Q-LdUID 

Ui  E  QQQ 

LlI  LlI  LlI  LlJ  ■ 


5-mmm 

<Z 

tOU.  U  U  Q 

—J 

M 

tooooa: 

<r 

UI  UJ  UI 

3> 

E  CL  CL  Q-  tO 

u>  j 

<r 

on  to  to  to  ui 

O  Q 

Z! 

to 

UCklCLO; 

■  M 

z 

►"•uiuuiz 

ocet 

U<  Ui 

o 

z  to  to  to© 

h  fr~ 

ui  y 

►-« 

D3DDH 

Ou 

£ 

Si 

o 

w-icuco^r  o 

39 
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PLACE  CURSOR  OH  THIS  LINE  AND  PRESS  PROPER  KEY  FOR  DESIRED  ACTION 
KEY  ACTION 

1  JOB  TERMINATION 

2  NEXT  CASE-NEU  A  AND  B 

3  NEXT  CASE-OLD  A  AND  B 
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THE  PRIOR  BETA  PARAMETERS  FOUOU 
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THE  MEAN  OF  RS  ■  *85277 
THE  MODE  OF  RS  ■  .86097 
THE  UAR  OF  RS  •  .00141 


n  m 


3  3  3 

u>  •  « 


o 


r* 
.  83 
Z  isl 
8  sf 

IA 


hi 


o 

u. 

V* 


*  s 

u  v> 

*  s 

Q 

S  2 

t- 

Ui 

c 

8  J, 


*-<  m 
t  r- 


l 


2  J5 

Ul 

a 

£  ■ 
8  S 


2  S  £ 

y  a  a 


Ui  Ul  U| 

sh 

it  £  £ 


f  ? 


S  • 


K  ft  dC 

3  3  3 

o  o  o 

J  J  J 

it  it  it 


o 

o 


it 

o 

8 


tx> 

r> 


je  20 


o 

c 


o 

U- 

tn 


«noMu 

M9A 

<v£ 

V'  V) 
00  00 
u>  •  • 

m 

to 


o  r-  •  • 
x  r- 
Ul  V>  v> 

a  at  a 

u  »  ru 
z  4o 

ss 


Ul  Ul  Ui 

ss  ss  ss 

I  I  I 


00 

o 

o 

o 

▼ 

H 

h 

1- 

f\J  » 

x 

x 

00 

m  v> 

Ul 

V 

ui 

Ui 

it  it  it 

XXX 
h  h  K 

S  8  K 

ft  ft  ft 

3  2  3 

j  j 

it  it  it 


S 

r- 

ai 


o 

x 


o 

o 


O  ^  O 


45 


PLACE  CURSOR  ON  THIS  LINE  AND  PRESS  PROPER  KEY  FOR  DESIRED  ACTION 
KEY  ACTION 

1  JO»  TERMINATION 

2  NEXT  CASE-NEW  A  AND  I 

3  NEXT  CASE-OLD  A  AND  9 
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THE  PRIOR  BETA  PARAMETERS  FOUOU 
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APPENDIX  A 

MANUAL  COMPUTATION  FOR  A 
SIMPLE  TEST  CASE 
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MANUAL  COMPUTATIONS  FOR  A 
SIMPLE  TEST  CASE 


The  benefits  of  BIGRAP  can  be  readily  demonstrated  by 
comparinq  the  ease  with  which  BIGRAP  is  used  with  the  cor¬ 
responding  alternative  manual  computations.  In  this  appendix 
the  manual  computations  required  to  perform  just  the  Method  of 
Moments  Bayesian  reliability  assessment  for  a  simple  system 
configuration  will  be  compared  to  utilization  of  BIGRAP. 

The  reader  is  reminded  that  BIGRAP  also  provides  the  Monte 
Carlo  option  for  which  manual  computations  are  virtually 
impossible . 
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Example : 


a.  Reliability  Block  Diagram: 


b.  Boolean  Expression: 


A*B* (C+D) 


(Note:  For  more  complex  systems  conversion  of 

Boolean  expression  to  algebraic  expression  is  a 
cumbersome  task  without  the  aid  of  a  computer 
program) . 

c.  Component  Descriptions: 


Name  Type , 


Successes , s. 


Data 


Failures,  f. 


A 

B 

C 

D 


R] 

R: 

R. 


R. 


19 

17 

35 


1 

3 


R  : 
s 


d.  Algebraic  Reliability  Expression  for  System  Reliability, 


Rs  2R1R2R3  rir2R3 


e.  Square  of  Algebraic  Expression: 

R?  =  4R2R2R2  -  4R?R2R2  +  R2R2R^ 
s  123  123  123 


f.  First  Moment  of  System  Reliability  Distribution  in 
Terms  of  Moments  of  Component  Reliability  Distributions: 


E(Rg)  =  2E(R1)E(R2)E(R3)  -  E(RL)E(R2)ECR3) 
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Second  Moment  of  System  Reliability  Distribution: 

E  (Rg)  =  4E  (R^ )  E  ( R^ )  E  ( R^ )  -  4E  (R^)  E  (R^)  E  (R^) 

+E  ( R^ )  E  ( R^ )  E  ( R^ ) 

h.  Find  Component  Prior  Beta  Parameters,  a'  and  b'. 

Which  Yield  Uniform  Prior  Distribution  for  System  Reliability: 


(1)  Substitute  E(R^)  =  .n  r^.  fj-T ^ 

1  3=1  a  +b  +3— 1 

for  each  expectation  appearing  in  first  and  second  moment 
expressions  above. 

(2)  Set  resulting  first  moment  expression  equal  to  1/2, 
the  value  of  the  first  moment  of  a  uniform  distribution  on  the 
interval  (0,1) . 

(3)  Set  resulting  second  moment  expression  equal  to 
1/3,  the  value  of  the  second  moment  of  the  uniform  distribution. 

(4)  Solve  the  resulting  two  equations  for  the  two 
unknowns,  a'  and  b'.  For  this  simple  system  the  two  equations 
to  be  solved  are: 


E(V 


2  r_2J _ } 

la '  +b '  > 


r  a'  ' 3  a 1 +1 
'•a'+b''  a'+b'+l 


1/2 


and 


E(Rs) 


,  r  a'  )3r  a'+l  ^ 3  ,  a*  ^  3  (  a'  +  l  ^  3  a '  +2 

la'+b,;  '■a'+b'+l'  ^a’+b’^  ^ 1  +b  •  +P  a  •  +b  •  +2 


+ 


a  ’  +b 


^)3( 


l±i  )  3r 

•tb'+l7  1 


a' +2  w  a* +3  ) 

a,+b,+2Ma,+b,+3/ 


1 

3 


To  solve  these  non-linear  equations  without  the  use  of  the 
computer  is  a  difficult  and  time  consuming  task.  Iterative 
procedures  such  as  Newton's  Method  must  be  employed. 
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The  solutions  for  this  example  are: 

a'  =  1.2382 
b'  =  0.3965 


The  values  are  used  as  the  prior  beta  parameters  for  all 
components  and  result  in  an  approximately  uniform  prior 
distribution  for  system  reliability. 

i.  Component  Posterior  Beta  Parameters: 

al  =  ®1  +  a'  =  19  +  1.2382  =  20.2382 

bl  =  fl  +  b'  =  1  +  0.3965  =  1.3965 

a2  =  s2  +  a'  =  17  +  1.2382  =  18.2382 

b2  =  f  +  b'  =  3  +  0.3965  =  3.3965 

a3  =  s3  +  a'  =  35  +  1.2382  =  36.2382  • 

b3  =  f  +  b'  =  5  +.0.3965  =  5.3965 

j.  Component  Posterior  Moments  Needed  in  First  and 
Second  System  Moment  Expressions: 


EU^)  = 


al+bl 


20.2382 

21.6347 


=  0.93545 


E(Rp  = 


al  al+1 

a^+b^  *  a^+b^+1 


=  0.87774 


E(R2) 


a2+b2 


=  0.84301 


e<R2> 


a2+b2 


a2+1 

a2+b2+1 


=  0.71651 
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0.87038 


E(R-)  =  — fr— 

3  a3  b3 


e,r3>  -  e<r3>- 


=  0.76022 


3  3  a~+2 
E(R3>  "  E<R3>*  a^+b^+2 


=  0.66620 


4  3  a->  +  3 

E(R  )  =  E(Rf).  fr- 
3  3  a2+b^+3 


=  0.58565 


k.  First  Posterior  Moment  (Mean)  of  System  Reliability: 

Substitute  the  results  in  j .  in  first  moment  expression 
provided  in  f. 

E (Rg)  =  2  (0.93545)  (0.84301)  (0.87038) 

-  (0.93545)  (0. 84301)  (0. 76022) 

=  0.77325 

l.  Second  Posterior  Moment  of  System  Reliability: 

Substitute  the  results  in  j .  in  the  second  moment 
expression  provided  in  q. 

E(R2)  =  4  (0.87774)  (0.71651)  (0.76022) 

-4  (0.87774)  (0.71651)  (0.66620) 

+  (0.87774)  (0.71651)  (0.58565) 

=  0.604851 


m.  Posterior  Variance  of  System  Reliability: 
Var (Rg )  =  E(R2)  -  E2(Rg) 

=  0.604851  -  (0 . 77325) 2 
=  0.006936 
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n.  System  Reliability  Posterior  Beta  Parameters,  a  and  b  : 

s  s 

E2(R  )  ( 1-E (R  )  ) 

Use  a  =  - - - - - 

Var(Rs) 

E ( R  )  ( 1-E ( R  )  )2 

,  _  S _  S 

s  Var(R  ) 

s 

Hence  ag  =  18.77365 
bg  =  5.50523 

o.  Beta  Density  Function  of  Posterior  Distribution  of 
System  Reliability: 


-  E<V 

+  E(Rs)-l 


f  B ( Rg|  as,bs)  =  f  ^ ( Rj  18.77365,5.50523) 

r(24.27888)R  17‘77365(1-R  )4* 50523 
_  _ _  s _  s 

r ( 18 . 77365)  r (5.50523) 


p.  Best  Estimate  of  System  Reliability: 

In  this  case  the  posterior  distribution  has  a  mode 
less  than  1.0  so  that  the  mode  is  the  best  estimate.  * 

a  -1 

Mode  of  R  =  — xr - a-  ~  0.79778  (Best  Estimate) 

s  a  +b  -2 
s  s 

q.  100  (1-Y  ) %  Lower  Confidence  Limit  on  System  Reliability: 
Find  R  such  that 

R 

/  f  (R  |  a  ,b  ) dR  =  Y 
J  s  s 9  s  s 

o 

(Note:  Here  aaain  manual  computations  are  difficult  requiring 

complex  interpolation  schemes  of  tabled  incomplete  beta  values 
or  computer  numerical  integration  algorithms) 

*See  page  15  for  explanation  (underlined  line) 
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APPENDIX  B 

THE  SIMPLE  TEST  CASE 
AS  RUN  WITH  BIGRAP 
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UITH  USER-ID  PG 
EOUIP/PORT  40/006 
COIWAND-  ETL.300 

COIWAND-  BEGIN,  RELY,  BAYE/CIARDO,  XX,  YYY 


ENTER  TOTAL  NUMBER  Of  DIFFERENT  COMPONENT  TYPES 
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THE  PRIOR  BETA  PARAMETERS  FOLLOW 
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The  numerical  differences  between  the  results  of 
the  computer  program  and  the  manual  computations  can 
be  attributed  to  the  fact  that  the  calculations  in  the 
computer  program  are  based  on  an  numerical  approxima¬ 
tion  technique. 


STOP  4 

602  CP  SECONDS  EXECUTION  TINE 
284  CP  SECONDS  COMPILATION  TIME 
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20,  ABSTRACT  (Vcmitnue  en  rereree  et da  ft  neceaeary  mtd.  Identify  by  block  number) 

BIGRAP  is  a  package  of  interactive  graphics  proqrans,  written 
for  use  on  the  TEKTRONIX  4014  Graphics  Terminal  connected  to  the 
ARRADCOM  CDC  6500/6600  computer  configuration.  This  nackane  con¬ 
sists  of  a  set  of  intricate  oroarams  that  allows  a  user  to  input 
component  success/failure  data  and  a  Boolean  expression  depictina 
system  reliability  locric  for  the  purpose  of  assessino  svstem 
reliability.  The  computer  converts  the  Ionic  expression  to  an  • 
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20.  Abstract  (continued) 

alqebraic  exnression  For  the  system  reliability  as  a  ^unction 
of  the  individual  component  reliabilities.  A  Bayesian 
statistical  alqorithm  is  them  employed  to  provide  the  user 
with  point  and  confidence  interval  estimates  of  system 
reliability.  In  addition  the  araohics  feature  o^  the 
package  displays  histoqrams  and  correspondinq  Beta  distri¬ 
butions  involved  in  the  analysis.  The  oroqrams  are  run 
sequentially  as  a  cataloaed  nrocedure  usina  the  BEGIN/REVET'.T 
utility.  This  reoort  is  intended  to  serve  as  a  quide  to 
potential  users  of  RTGRAP  for  operation  and  interpretation 
of  results. 


